An efficient approach of controlling traffic congestion in scale-free networks 
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We propose and study a model of traffic in communication networks. The underlying network has a struc- 
ture that is tunable between a scale-free growing network with preferential attachments and a random growing 
network. To model realistic situations where different nodes in a network may have different capabilities, the 
message or packet creation and delivering rates at a node are assumed to depend on the degree of the node. 
Noting that congestions are more likely to take place at the nodes with high degrees in networks with scale-free 
character, an efficient approach of selectively enhancing the message-processing capability of a small fraction 
(e.g. 3%) of the nodes is shown to perform just as good as enhancing the capability of all nodes. The interplay 
between the creation rate and the delivering rate in determining non-congested or congested traffic in a network 
is studied more numerically and analytically. 

PACS numbers: 89.75.Hc, 05.70.Hc 
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INTRODUCTION 

Operations in the internet such as browsing webpages in the 
World Wide Web (WWW), sending e-mails, transferring files 
via ftp, searching for information, and electronic shopping, 
etc. have become part of daily life for many people. These 
activities have opened up exciting opportunities for sharing 
information, economic transformation, and other activities on 
a global scale |Q1 01 . The internet, however, is not perfect. 
For example, intermittent congestion in the internet, similar to 
traffic congestion in highway systems, has been observed [3]. 
Similar phenomena can be also of relevance in other commu- 
nication networks, such as the transportation network in air- 
lines and the postal service network. A key problem in com- 
munication networks is, therefore, to understand how one can 
control congestion and maintain a normal and efficient func- 
tioning of the networks. 

Several models of communications in a computer network 
have been extensively studied flBllBSSGSOilEi 
. 14. 15. [lrjll . In these models, the information proces- 
sors are routers. Their function is to route the data packets 
to their destinations. In a computer network, a node may be 
a host or a router. A host can create messages or data pack- 
ets for targeted destinations and receive packets from other 
hosts. A router finds the shortest path between the origin of 
the message and the destination of each packet and forwards 
the packet one step closer to the destination along the shortest 
path in each time step. The shortest path is the path with the 
smallest number of links. Previous studies have mostly been 
focused on three different computer network models: (i) the 
nodes at the edge of the network are hosts and the inner nodes 
are routers |U, (ii) all the nodes are both hosts and routers 
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, and (iii) some of the nodes are hosts and the rest 
are routers ifTH [l5l [lrjll . However, these models were studied 
with the underlying networks being a two-dimensional lattice 
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net shows a het erog eneous structure with a scale-free degree 
distribution 18], a more realistic network model for com- 
munications should be heterogeneous. Besides the internet 
and its related networks such as WWW Bill and email 

networks J22I1 . many other networks also show the scale-free 
behavior. These networks include, for example, the telephone 
network 12311 . the biological network in proteins B24I1 . and the 
networks of sexual contacts [25, 26j]. Some networks, such 
as the collaboration network among scientists ll27ll . show a 
mixed feature of scale-free and exponential distributions. In 
fact, the study of the science of complex networks has become 
an important interdisciplinary area of research. The problem 
of efficiency in delivering messages or data packets in com- 
munication networks has been addressed recently by Arenas 
et alJM HI |3(1, Moreno et al. JH 0,0], and Zhao et 
al. 113411 . Arenas et al. focused on finding the optimal network 
topologies for searches in complex networks, while Moreno et 
al. studied the dependence of the jamming transitions on rout- 
ing strategies. A common feature in previous studies is that 
the creation and delivering rates of packages do not change 
from node to node. As the nodes in a complex network could 
have very different properties, e.g., degrees, a more realistic 
assumption is that the package creation rate and delivering 
rate at a node become degree-dependent. In the internet, an 
important site has more users and hence a larger message or 
package creation and delivering rates. A recent study by Zhao 
et al. IB4I1 considered the case of non-uniform package deliv- 
ering rates, but the creation rate was taken to be a constant. 

In the present work, we study traffic in networks with non- 
uniform package creation and delivering rates. An important 
quantity in communication networks is the critical package 
creation rate that signifies a transition from a non-congested 
or free flow regime to a congested regime. Below the critical 
rate, a non-congested steady state is reached after the tran- 
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sient in which the data packets created can be efficiently han- 
dled by the nodes. Above the critical rate, a congested phase 
is reached where the number of packets accumulated in the 
system increases with time. The value of critical rate thus 
measures the capacity of efficient communication inside the 
network. Here, we study the critical rate in networks where 
the package creation and delivering rates are node-dependent. 
In particular, we present an efficient approach to enhance the 
capacity of communications in scale-free networks. 

The paper is organized as follows. Section II defines the 
model with node-dependent package creation and delivering 
rates. In Section EH, we present results of numerical simula- 
tions and study the interplay between the critical and deliv- 
ering rates in determining non-congested or congested traffic 
in a network. Section IV explains the observed features in 
the numerical results analytically. We summarize the paper in 
Sec. V. 



MODEL 

The nodes in a complex network such as the internet may 
represent very different entities. For example, some nodes 
may just be individuals and other may represent big com- 
panies or universities. Obviously, different nodes will have 
different rates of creating messages. The nodes, depending 
on their connectivity to other nodes and perhaps hardware, 
also have different rates of delivering messages. Here, we 
present a more realistic model of communication in complex 
networks that includes node-dependent creation and deliver- 
ing rates. Our model is a modification on several previous 
models JSISIMIIIIIlSB.We assume that for a node 
i with degree ki, the message creation rate Xki is proportional 
to its degree, with A being a constant. For message deliver- 
ing, each node should handle at least one packet or message 
in each time step. Therefore, we assume a delivering rate of 
1 + (3ki for a node with degree ki, with f3 > being a pa- 
rameter of the model. Our model thus represents the realistic 
situation that a busy node with larger ki has higher rates of 
generating and delivering messages. 

For the underlying network, we use a model in which the 
exponent of the degree distribution can be tuned. A scale- 
free network with P(k) ~ fc -7 can be constructed by incor- 
porating preferential attachments in a network-growing pro- 
cess 11351 13611 . The Barabasi and Albert model H36I1 assumes 
the probability for a node i to attract a link from a newly 
added node to be ~ fcj. The model gives an exponent 
7 = 3 ll36ll for the degree distribution. On the other hand, a 
random growing network can be constructed by assuming a 
node-independent n^. Many networks show characters that 
are somewhat intermediate of scale-free and random. For 
them, the degree distribution shows a mixed feature of the 
two characters 112711 . This implies that the probability of 
attracting a new link should contain both preferential and ran- 
dom features. One of the present authors proposed a hybrid 
model [37] in which ~ (1 — p)ki + p, where < p < 1 is 



a parameter representing the probability that a newly added 
node establishes its new links by random attachments and 
(1 — p) is the probability that new links are established by 
preferential attachments. The degree distribution was shown 
to be Q P(k) ~ [k + p/(l - p)]" 7 ^ with an exponent 
7(p) = 3 + p/\m{l — p)], where m is the number of new 
links per node. The p = limit reduces to the P(k) ~ fc~ 3 
behavior and the p — * 1 limit gives the random growing net- 
work behavior of P{k) ~ e~ fc / m . Here, we use this model as 
our underlying network for studying communications in net- 
works. 

Once the network of a certain value of p is constructed, the 
dynamics of creating and delivering messages is implemented 
as follows. Each node plays the dual role of a host and a 
router, with its creation and delivering rates assigned accord- 
ing to its degree. Details of the dynamics are listed as follows. 

(1) At each time step, a node i has a probability Xki of cre- 
ating a new message or packet with a randomly chosen desti- 
nation. If the node has some messages waiting to be sent, the 
newly created message will be placed at the end of the queue. 
The queuing messages may be created at some previous time 
steps or received from nearby nodes as messages are being 
sent along their paths to the destinations. 

(2) Once a packet is created with a chosen destination, the 
node (router) will identify the shortest path towards the desti- 
nation. If there exist several shortest paths to the destination, 
the path is chosen in such a way that the package is sent to a 
node that has the instantaneous shortest queue. 

(3) At each time step, a node i has the ability to forward 
( 1 + j3ki ) packets in the queue at the node on a first-in-first-out 
basis to its neighbors which are along the path to the destina- 
tions. Noting that f3ki may be an integer plus a fractional part, 
the fractional part is implemented as the probability of deliv- 
ering additional packets in a time step. 

(4) Messages arriving at a node are queued up for further de- 
livering. When a message arrives at its destination, it is re- 
moved from the system. 

The steps are carried out for every node at the same time. If 
Xki is replaced by A and (3ki is replaced by the integral part 
int[(3ki], the above algorithm will be equivalent to that of Ref. 
lT34h . Here, the fractional part of (3ki and Xki are implemented 
in a probabilistic way. Furthermore, if we take (3 = 0, the 
above algorithm will be equivalent to that of Refs j4 nan 
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The parameters A and (3 thus control the number of mes- 
sages or packages in the system. A small (high) value of A cor- 
responds to fewer (more) packets. For simplicity, each packet 
is labeled by two pieces of information: the time of creation 
and its destination. Qualitatively, the total number of pack- 



ets created at each time step is J^iLi ~ $2»=i Mty 
2mXN for a growing network of m newly added links per 
node and a total of N nodes. At the same time, the maximum 
number of packets processed by the nodes is ^^ =1 (l+/3fcj) ~ 
(1 + 2mf3)N. When the number of new packets added to the 
system equals the number of packages removed upon arrival at 
each time step, the network runs in the range of Little law 
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and there is no congestion. For scale-free networks, the struc- 
ture is heterogeneous in the sense that some nodes have many 
more links. As messages are sent via the shortest paths, they 
are likely to pass through the nodes with more links. If every 
node has the same delivering rate, these nodes will have more 
accumulated messages and a high chance of jamming. In con- 
trast, random networks do not have hubs with high degrees 
and thus the structure is relatively "homogeneous". There- 
fore, congestion is easier to occur in scale-free networks than 
random networks. As most of the real-life networks are scale- 
free, the control of congestion in these networks is of critical 
importance. For example, an intuitive but not so efficient ap- 
proach is to increase the value of (3 for the nodes. We will 
refer to this approach of having identical values of (3 for all 
nodes as the normal approach. 

The number of nodes in a communication network is typ- 
ically very large and there is no central organizer to man- 
age the development of the whole network. It is, therefore, 
very difficult to have a realistic mechanism to increase (3 for 
all nodes at the same time. Realistically, larger companies 
and academic institutions could increase their local value of (3 
more readily. Noting that congestion is more likely to occur 
at the nodes with many links in scale-free networks, one may 
significantly reduce congestion by selectively increasing the 
delivering rates of the nodes with high degrees. Therefore, we 
suggest that a possible cost-effective way to control conges- 
tion is to ask the nodes with larger links to increase their value 
of (3. Here, we study a model in which a fraction / of nodes 
with high degrees are assigned a finite value of (3 > 0, and 
the rest are assigned (3 — 0. This models the higher message- 
processing capability of the hubs in a network. We refer to this 
model as the efficient approach. In this model, the maximum 
number of packets processed in a time step is (iV+$^ e / (3ki), 
where the sum is over the nodes with finite (3. Results of nu- 
merical simulations show that the efficient approach performs 
comparably with the normal approach. In the following sec- 
tions, we compare results of the two approaches and explain 
the results analytically. 



NUMERICAL RESULTS 

The network is constructed as described in Sec. II and in 
1000, m = 3 and different values of p. 



Ref.[37], with N 
The dynamics of package delivering is then implemented on 
the network. We first consider the normal approach in which 
all the nodes have the same value of (3. Intuitively, a larger 
(3 can assure free traffic flow for a larger creation rate signi- 
fied by a larger value of A. Here, we fix A = 0.01 and take 
(3 = 0, 0.05, and 0.1 to illustrate the effects. For (3 = 0, every 
node has a creation rate Afc; depending on fej, but the deliver- 
ing rate of forwarding at most one message per time step ap- 
plies to all nodes. To understand how congestion occurs, we 
calculate the average number of packets (n(k)) on the nodes 
with a given number of links k. This quantity serves to show 
where are the longest queues. Numerical simulations show 




FIG. 1: The average number of packets (n(k)) as a function of 
the number of links k in networks characterized by m = 3 and 
N = 1000 for (a) random growing networks (p = 1) and (b) scale- 
free networks (p = 0). The parameter characterizing the message 
creation rate is A = 0.01. Results are obtained after t = 500 time 
steps and averaging over 100 different realizations for a given set of 
parameters. Different symbols label different packet-processing ca- 
pabilities: ft = (circles), (3 = 0.05 (squares), and (3 — 0.1 (stars). 



that they are the nodes with more links, as shown in FigQ](cir- 
cles) for systems after t = 500 time steps. The results are ob- 
tained by averaging over 100 different realizations for a given 
set of parameters. Fig[TJa) shows the results for random grow- 
ing networks (p = 1) and FigJTtb) shows the results for scale- 
free network (p = 0). It is clear that in both cases, the nodes 
with large number of links are more likely to be congested. 
Comparing the results in FigOJa) with (b), one sees that the 
accumulation of packets in scale-free networks is much pro- 
nounced than that in random networks. The results indicate 
that congestion is much easier to occur in scale-free networks. 
For (3 = 0.05 (squares in FigJTJ, the accumulation of packets 
is greatly suppressed in both limits of the underlying network. 
In particular, congestion almost disappeared in the case of ran- 
dom growing networks. For even higher message-processing 
capability (3 = 0.1 (stars in FigJTJ, congestion disappeared in 
both the random and scale-free networks. These results also 
indicate that there exists a critical value (3 C for a given A so 
that congestion occurs for (3 < (3 C . We will study the depen- 
dence of (3 C on A in networks of different underlying structures 
characterized by p. 

Another point that is worth noticing is that the nodes in a 
range of small to intermediate degrees (see inset in Fig[T|) in 
the normal approach usually carry fewer messages than the 
uniform delivering capacity. The result implies that it is un- 
necessary for these nodes to have a higher delivering capa- 
bility characterized by a finite (3. It leads us to consider the 
efficient approach. In the scale-free (p = 0) limit as shown in 
Figffib), t ne degree distribution is a power law and the nodes 
that carry k > 20 account for only 3% of all the N = 1000 
nodes. To illustrate the idea of the efficient approach, we take 
/ = 3%, i.e., we assign a non-vanishing (3 only to nodes 
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FIG. 2: (n(k)) as a function of k for scale-free networks (p = 0) 
with TV = 1000 and m — 3 for the efficient approach (circles) with 
/ = 3% and the normal approach (stars). The parameter characteriz- 
ing the message creation rate is A = 0.01. Results are obtained after 
t = 500 time steps and by averaging over 100 different realizations 
for a given set of parameters. Two values of /3 are used: (a) j3 — 0.05 
and (b)/3 = 0.1. 



FIG. 3: (a) The average number of messages per node and (b) the av- 
erage number of messages among the top 3% nodes with the highest 
degrees as a function of time in a scale-free network (p = 0) for three 
different values of /3. The lines from top to bottom refer to (3 = 0.05, 
0.059, and 0.07, respectively. The other parameters are A = 0.01, 
/ = 3%, N = 1000, and m = 3. 



that have k > 20. Figure [2] compares results of the normal 
approach (stars) and the efficient approach (circles) for two 
values of (3. Obviously the difference in (n(k)) between the 
two approaches is small in scale-free networks. As the effi- 
cient approach does not require all nodes to be equipped with 
the same capability, it represents a more practical and cost- 
effective way to avoid jamming. In contrast, we note that the 
efficient approach does not work so well in random networks 
(p = 1). It is because the degree distribution is narrower com- 
pared with the p = case. Therefore, the queues are more 
evenly distributed among the nodes and congestion is not re- 
stricted to the nodes among the highest degrees. It is thus 
necessary to assign a finite / to a larger fraction of nodes to 
avoid congestion, and the two approaches become similar. 

We define (rii (t)) to be the average number of messages per 
node. In the congested regime, (m(t)) increases with time 
t. In the non-congested regime, (n 1 (t)} fluctuates around a 
constant. The slope of < (t) > after the transient can thus 
be used to determine /3 C OOi 13211 . For a given A, the slope 
gradually decreases as (3 increases. The value of (3 that the 
slope becomes zero gives j3 c . For j3 > (3 C , the slope remains 
zero. Figure|3la) shows typical results with A = 0.01 for three 
different values of (3 within the efficient approach (/ = 3%). 
For (3 = 0.05, (m(t)) increases with time without bound. 
The critical value is found to be (3 C = 0.059 where the slope 
vanishes. For (3 — 0.7 > j3 c , the slope remains zero. As 
congestion mainly occurs at the nodes with large degrees, the 
number of messages (n 2 (t)) averaged over the 3% of nodes 
should also show a similar behavior with time. It is indeed the 
case (see FigJ5Jb)). 

Next, we study the dependence of f3 c on the creation rate 
characterized by A in scale-free and random growing net- 
works. FigUa) shows the results in the scale-free limit 



(p = 0) for both the normal (circles) and efficient (stars) 
approaches. For small A, (3 C vanishes as the default deliver- 
ing rate of one message per time step is already sufficient to 
handle the small message creation rate. For the range of A 
shown in the figure, (3 C increases linearly with A and the two 
approaches give similar results. This again shows that the ef- 
ficient approach performs as good as adjusting (3 across the 
whole network. It should be noted that for larger values of A 
(beyond the range shown here), assigning a finite (3 to only 
the top 3% of nodes may not be sufficient to avoid conges- 
tion. For random growing networks (p = 1), we show /3 C (A) 
in FigHfb) only for the normal approach, as the efficient ap- 
proach becomes similar to the normal approach. Qualitatively, 
/3 C (A) shows a similar behavior to that in scale-free networks. 
Quantitatively, (3 C — for a larger range of A in random 
networks and the slope of the linear dependence in f3 c (X) is 
smaller. It is because the nodes in random networks are more 
"homogeneous" and a queue will not emerge at the hubs for 
small A as in the case of scale-free networks. The function 
/3 C (A) in FigHfa) also divides the (3-X space into two regions. 
The region above the line represents a non-congested or free 
flow regime and that below the line represents a congested 
regime. Thus for given A, one can go from a congested to a 
non-congested regime by increasing [3. Similarly, for a given 
j3, one can go from a non-congested regime to a congested 
regime by increasing A. Although we only present results for 
networks with N = 1000 nodes, we have checked that the lin- 
ear dependence of (3 C on A also holds for networks with larger 
N. 
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FIG. 4: /3 C as a function of A for (a) scale-free networks (p = 0) of 
N — 1000 and m = 3 within the normal (circle) and efficient (stars) 
approaches (/ = 3%); and (b) random growing networks (p = 1). 



THEORETICAL EXPLANATION 

In this section, we aim at explaining the behavior of /3 C (A) 
quantitatively. As discussed, (n\(t)) behaves differently for 
(3 < j3 c and/3 > (3 C . In the non-congested or free flow regime, 
the steady state satisfies the Little's law 13811 . which states that 
the number of delivered messages is balanced by the number 
of newly created messages. This suggests a way to estimate 
/3 C (A) for a given underlying network characterized by the pa- 
rameter p. 

Consider the node with the highest degree k max (p) where 
messages are most likely to accumulate. At the critical value 
(3 C , the node can handle 1 + (3 c k max (p) packets per time step, 
while the number is smaller for other nodes. For scale-free 
networks (p = 0), k max w mTV 1 ^ 7-1 ) We note that 

the packets at a node i originate from two different sources: 
those created at node i and those passing by node i. The cre- 
ation rate Xki is linear in fcj. The packets passing by are more 
likely to go through the nodes with higher degrees, and hence 
the number of packets passing by a node will be some nonlin- 
ear function of its degree. With these considerations, we ap- 
proximate the average number of packets at the critical value 
f3 c at some node i as a(ki,p)(l + /3 c k max (p))ki/k max (p), 
where < a(ki,p) < 1 and a(k max ,p) = 1 is a nonlin- 
ear decreasing function of fej that reflects the contribution of 
messages passing by the node. For the case of p = 0, noting 
that there are only 1 + (3 c k max (fi) packets at the nodes with 
kmax, the average number of packets at the nodes with small 
and intermediate degrees will be less than one. This implies 
that there is not enough packets for the parameter (3 to take 
effect at these nodes. Therefore, we expect the expression 
a {h , p) (1 + p c k max (p))ki/k max (p) to be a good approxima- 
tion for both the normal and efficient approaches in scale-free 
networks. On the other hand, there are 2mXN newly created 
packets in each time step. Let h(p) be the diameter of the 
network which measures the average number of nodes that a 
packet passes through on its way to its destination, including 



the destination itself. If the system is in the non-congested 
regime, there are a total of h(p)2m\N messages in the sys- 
tem. To avoid a queue at any node and hence congestion, all 
the messages should be handled by the nodes in a time step. 
Thus, we have 

N 

h(p)2m\N = y^a(kj,p)(l + (3 c k max (p))h/ k max (p). 
i=i 

(1) 

Writing Ei=i"(*i.P)(l + max (p))h/k (p) = 
a i(p) YsiLii 1 + Pckmax{p))ki/k rnax (p), we then have 



ai{p) K- n t a x (pY 

From Eq.©, it follows that (i) (3 C > only when A is suffi- 
ciently large, (ii) for a given structure of the network (fixed p), 
(3 C increases linearly with A, and (iii) the slope of /3 C (A) de- 
pends on the underlying network structure characterized by p. 
All these features agree with those observed in the numerical 
results (see Figi4j». 

Equation (O can be applied to estimate (3 C , if we know 
h(p), k max (p), and ai(p). The diameter h(p) can be cal- 
culated using the method in Ref. 113311 . Figure [5j a) shows 
h(p) over the whole range of p. It increases only slightly as 
p increases. On the other hand, k max (p) drops sensitively 
with p as shown in Fig. |3b). Since a\(p) depends only 
on p, it can be determined by using numerical results of (3 C 
for given A. For example, /3 C (A = 0.01) = 0.059 in the 
scale-free (p = 0) limit. Together with h(Q) — 3.32 and 
k max (0) = 85 (see Fig©, Eq.© gives ai (0) w 0.4522 and 
the slope of the line C (X) is h(p) / a\(p) = 7.34. Similarly, 
/3 C (A = 0.012) = 0.027 in the random network limit (p = 1). 
Together with h(l) = 3.82 and k max (l) = 25 (see Fig©, 
Eq.© gives ai(l) ~ 0.6842 and the slope of the line /3 C (A) 
to be 5.58. These values are in reasonable agreement with 
the slopes in the plots in Figj4] Equation (O also shows that 
(3 C = for A < X min = a 1 (p)/[h(p)k max (p)}. Using the 
extracted values of the parameters, we get A m i„ = 0.0016 
forp = and \ m ; m — 0.0072 for p = 1. These values are 
consistent with the results in FigH] 

CONCLUSIONS 

In network communications, a simple way to control net- 
work traffic is to limit the length of the queues M4011 . e.g. 
by source quenching, random dropping, fair queueing, etc. 
This will, however, increase the average delivering time. 
As many real-life networks are heterogeneous networks and 
many shortest paths between any two nodes pass through the 
nodes with high degrees, it will be these nodes that control 
the network traffic. Thus we study the strategy of enhancing 
the message delivering capability selectively at the nodes with 
high degrees. We found that the strategy works well in net- 
works with scale-free character and it is a highly cost-effective 
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FIG. 5: (a) The diameter h(p) and (b) the highest degree k max (p) in 
a network as a function of the underlying network structure charac- 
terized by the parameter p. The networks have N = 1000 nodes and 
m — 3. 



way to avoid network congestion. This idea is in line with the 
recent results in Ref. 114 lL 14211 . 

The major difference in network congestion in a scale-free 
network and a random growing network is that the scale-free 
network has hubs, i.e., nodes that are connected to many other 
nodes. The degree distribution in a scale-free network fol- 
lows a power law for large networks. In a random network, 
the degree distribution is relatively narrower and the degrees 
of the nodes do not differ by much. For identical message 
creation rate and delivering rate at the nodes, it is then ex- 
pected that congestion will take place mostly at the nodes of 
high degrees in a scale-free network. For a random network, 
congestion may take place at more places across the network. 
Strategically enhancing the message-processing capability at 
the high-degree nodes in a scale-free network as in the ef- 
ficient approach studied in the present work will greatly en- 
hance network traffic. This strategy also makes good use of 
the power-law degree distribution in that it is sufficient to al- 
locate resources to enhance the capability of a small fraction 
of nodes with high degrees in a network in order to avoid traf- 
fic congestion. If we carry out the same strategy to a random 
growing network, a much larger fraction of nodes will be in- 
volved and hence the cost-effectiveness will be lowered. 

In summary, we have constructed and studied a model of 
communications in complex networks. We use a network 
model that can be tuned from the scale-free preferential grow- 
ing network limit to the random growing network limit. Our 
model assumes a message creation rate Xki that depends on 
the degree of a node. Each node also has a message delivering 
rates of 1 + (3ki. The model thus represents a step towards 
a more realistic modelling of traffic congestion in communi- 
cation networks in that it incorporates the different capacities 
of the nodes in creating and handling messages. In particular, 
we studied an efficient approach that increases the communi- 
cation capacity in scale-free networks. Numerical results in- 
dicate that our efficient approach of selectively enhancing the 



delivering rate in a small fraction of nodes performs as good as 
enhancing the capability of all the nodes in the network. Con- 
sidering the cost of enhancing the delivering rate at a node, the 
present scheme will be highly cost-effective. We also studied 
the dependence of the critical value of 0, which characterizes 
the message delivering rate, on the parameter characterizing 
the message creation rate A. The function /3 C (A) divides the 
j3-\ space into two regions of physically different characters: 
non-congested or free flow regime and congested regime. An- 
alytically, we derive an expression of f3 c (\) based on the idea 
that all the messages in the system should be handled by the 
nodes in the non-congested regime. The analytic expression 
captures all the features observed in numerical results. 
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